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Summary.  The  equilibrium  Joint  probability  distribution  of  queue  lengths  Is 
obtained  for  a  broad  class  of  Jobshop-Ilke  "networks  of  waiting  lines,"  where 
the  mean  arrival  rate  of  customers  depends  almost  arbitrarily  upon  the  number 
already  present,  and  the  mean  service  rate  at  each  service  center  depends 
almost  arbitrarily  upon  the  queue  length  there.  This  extension  of  the 
author's  earlier  work  Is  motivated  by  the  observation  that  real  production 
systems  are  usually  subject  to  Influences  which  make  for  Increased 
stability  by  tending,  as  the  amount  of  work- In-process  grows,  to  reduce 
the  rate  at  which  new  work  la  Injected  or  to  Increase  the  rate  at  which 
processing  takes  place. 


1.  Introduction 

This  paper  provides  the  equilibrium  Joint  probability  distribution  of 
queue  lengths  for  a  btoad  class  of  queueing-theoretical  models  representing 
multipurpose  production  systems  composed  of  special-purpose  service  centers 
(e .g ,  manufacturing  Jobshops) .  In  the  systems  modelled ,  customers  arrive 
from  time  to  time,  each  with  a  routing,  which  Is  an  ordered  list  of  service 
centers  where  the  customer  must  be  served.  An  arriving  customer  Joins  the 
queue  at  the  first  center  on  his  routing,  remains  there  until  his  service 
there  Is  completed,  then  goes  to  the  second  center  on  his  routing  (If  any) 
and  remains  there  until  served,  and  so  forth,  until  all  his  required  services 
are  completed,  at  which  time  he  leaves  the  system. 


1.  This  work  was  supported  by  the  Office  of  Naval  Research  under  Task 
047-003  and  by  the  Western  Management  Science  Institute  Ford  Foundation 
grant.  Reproduction  In  whole  or  In  part  Is  permitted  for  any  purpose 
of  the  United  States  Government. 

2.  To  be  presented  to  the  Joint  International  Meeting  of  the  Institute  of 
Management  Sciences  and  the  Operations  Research  Society  of  Japan,  Tokyo, 
August  21-24,  1963. 
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The  customer  arrivals  are  modelled  as  a  generalized  Poisson  process, 
whose  parameter  (mean  arrival  rate)  varies  almost  arbitrarily  with  the 
total  number  of  customers  already  in  the  system.  Service  completions  at 
each  center  are  also  modelled  as  generalized  Poisson  processes,  the  parameter 
(mean  service  rate)  at  each  center  varying  almost  arbitrarily  with  the  queue 
length  there  (number  of  customers  in  the  queue ,  waiting  and  being  served) . 

The  generation  of  routings  Is  modelled  as  a  random  walk  among  the  names 
of  the  service  centers ,  which  with  probability  one  reaches  a  terminal  state 
marking  the  ends  of  routings .  The  mechanism  for  generating  routings  Is 
sufficiently  general  to  Include,  for  Instance,  cases  where  every  conceivable 
routing  may  occur,  the  possibility  that  all  customers  require  services  at 
all  centers  In  the  same  order,  and  Instances  In  which  each  customer  requires 
service  at  just  one  of  the  centers. 

The  mathematical  model  does  not  explicitly  reflect  any  particular  rule 
for  deciding  which,  among  several  customers  at  a  service  center.  Is  to  be 
served  "next."  The  model  and  its  analysis  are,  In  fact,  consistent  with 
all  rules  which  do  not  depend  (directly  or  indirectly)  upon  the  future 
routing  or  service  time  requirements  of  any  customer.  Perhaps  the  most 
commonly  used  rule  In  this  category,  and  the  one  whose  operation  is  most 
easily  visualized,  is  the  first-come,  first-serve  rule. 

Results .  The  basic  result  of  this  paper  is  the  equilibrium  joint 
probability  distribution  of  the  queue  lengths  at  the  several  centers  in  the 
jobshop-like  queueing  system  just  outlined.  This  result  is  stated  as  a 
theorem  in  Section  4,  following  the  detailed  formulation  of  the  basic  model 
and  preliminary  analysis  in  Sections  2  and  3.  In  Section  5,  the  theorem  of 
Section  4  is  generalized  to  cover  cases  where  the  immediate  injection  of  a 
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new  customer  Is  triggered  when  the  total  number  of  customers  falls  below 
a  specified  limit,  or  where  a  service  Is  deleted  If  a  queue  grows  beyond 
a  specified  maximum  length.  Section  5  includes  remarks  about  the  theorem 
and  points  to  some  of  Its  implications. 

Section  6  is  devoted  to  special  cases.  Systems  with  constant  customer- 
arrival  rate  are  considered;  and  It  turns  out  that  the  main  theorem  of  Section 
5  reduces  in  this  case  to  a  statement  that  the  equilibrium  joint  probability 
distribution  of  queue  lengths  Is  identical  with  what  would  be  obtained  by 
pretending  that  each  Individual  service  center  is  a  separate  queueing 
system  Independent  of  the  others.  The  constant  arrival-rate  case  Is  further 
specialized  to  the  Instance  where  each  center  is  an  "ordinary"  multi-channel 

3 

service  system.  Next,  to  clarify  the  effect  of  varying  arrival  rates, 
attention  Is  directed  to  a  special  case  of  the  general  model,  in  which  all 
service  rates  are  the  same.  Finally,  systems  are  considered  In  which  the 
total  number  of  customers  present  Is  held  fixed. 

Significance.  The  generalized  customer-arrival  process  of  the  present 
paper  allows  for  the  representation  of  systems  whose  "potential  customers" 
are  generated  by  a  stable  Poisson  process,  but  where  the  probability  that  a 
potential  customer  will  actually  enter  the  system  for  service  depends  upon 
how  many  customers  are  already  there.  In  addition  to  the  representation  of 
"ordinary"  multiserver  centers,  the  generalized  service-completion  process 
also  permits  the  treatment  of  systems  at  whose  centers  the  number  of  servers 
actually  functioning  depends  almost  arbitrarily  upon  the  queue  length  there, 
as  well  as  plausible  rough  approximations  to  such  phenomena  as  "speedups"  and 
"overtime  work."  The  further  generalizations  (Section  5)  of  the  customer- 
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Which  was  analyzed  previously  in  the  author’s  paper,  "Networks  of  Waiting 
Lines,"  Operations  Research,  vol.  5,  no.  4  (Aug.,  1S57) . 
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arrival  and  service-completion  processes  make  possible  crude  approximations 
to  such  phenomena  aa  the  activation  of  ’’standby  Jobs”  when  a  shopload  falls 
off,  and  ’’subcontracting”  when  a  service  center  becomes  overloaded. 

2.  Notations  and  Restrictions 

The  first  business  of  this  section  is  to  specify  and  provide  notations 

for  the  parameters  of  the  generic  Jobshop-like  queueing  system.  As  this 

is  done,  these  parameters  are  interpreted  and  restrictions  are  stated.  The 

model  is  mathematically  formalized  in  Section  3. 

States  of  the  system.  The  number  of  service  centers  in  the  system  will 

be  denoted  by  N  ,  a  positive  Integer;  and  the  centers  will  be  referred  to 

as  Center  1,  Center  2,  etc.  The  states  of  the  system  are  N-dimensional 

vectors  with  non-negative  integer  components,  the  n-th  component  representing 

the  queue  length  at  Center  n.  Such  vectors  will  be  called  state  vectors. 

If  k  ss  (kj^,k2,  .  .  .  ,  kj^)  is  a  state  vector,  then  S(k)  =  k^^  +  k^  +  ,  .  •  +  » 

the  total  number  of  customers  in  the  system  when  its  state  is  k. 

4 

Customer-arrival  and  service-completion  processes.  The  cust<M&er-arrlval 
process  is  specified  by  a  set  of  parameters, 

L  =  {X(K)  I  K€  [0,co))  . 

The  interpretation  is  that  if  the  system  is  in  state  k  at  time  t  ,  then 
the  probability  that  a  customer  will  arrive  between  then  and  time  t+h 
is  h  \(S(k))  -f  o(h)  .  Similarly,  the  service-completion  processes  are 
specified  by  a  set  of  parameters, 

M  =  {u(n,10  1  n  €  [0,  ”))  ; 

With  the  interpretation  that  if  k  customers  are  present  at  Center  n 


4.  If  a  and  b  are  Integers  with  a  <  b.,  then  [a,oo)  denotes  the 
set  of  Integers  not  smaller  than  a  ,  and  [a,b]  denotes  the  set  of 
Integers  between  a  and  b  inclusive.  The  expression  o(h)  denotes 
unspecified  quantities  such  that  as  h  tends  to  the  limit  zero,  so 
does  o(h)/h  . 
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at  time  t  ,  then  the  probability  that  the  service  of  one  of  these 
customers  will  be  completed  by  time  t+b  is  h  (j,(n,k^)  +  o(h)  .  The 
probability  of  two  or  more  events  (arrivals  or  completions)  in  a  time 
Interval  of  length  h  is  o(h)  . 

Prom  their  interpretations,  the  X(K)  and  the  p,(n,k)  must  be 
non-negative,  and  also  each  p.(n,0)  =  0  .  However,  it  is  convenient,  and 
seems  to  result  in  no  loss  of  interesting  generality,  to  make  somewhat 
stronger  assumptions,  as  follows: 

ASSUMPTION  (2.1).  Either  the  > (K)  are  all  positive;  or  for  some 
non-negative  ,  X(K)  >0  «  K  <  \(K)  b  0  W  K  >  ^ 

ASSUMPTION  (2.2),  Tte  pi(n,k)  are  all  positive,  except  that 
each  p.(n,0)  s  0  . 

(These  assumptions  may  seem  to  exclude  the  possibly  Interesting  case 
where  a  service  center  is  "shut  down”  if  the  queue  length  there  falls  to 
some  specified  lower  limit.  This  case  can,  however,  easily  be  brought  into 
this  paper's  analytical  framework,  by  redefining  the  n-th  component  of  the 
state  vector  as  the  number  of  customers  at  Center  n  in  excess  of  the  lower 
limit  there,  and  appropriately  translating  the  arguments  of  \(K)  and 
p,(n,k)  ,) 

Generation  of  routings.  The  routing-generation  process  is  specified 
by  a  set  of  parameters, 

R  =  {r(m,n)  (  m  €[0,N]  ;  nG  [1,N}-1])  . 

The  interpretation  is  that,  for  m  and  n  €[l,N]  :  (i)  the  probability  is 
r(0,n)  that  Center  n  will  be  first  on  a  routing;  (ii)  the  probability  is 
r(m,Nfl)  that  if  Center  m  is  i-th  on  a  routing,  then  this  i-th  element  is 
the  last  one;  (ill)  the  probability  is  r(m,n)  that  if  Center  m  is  i-th 
on  a  routing,  then  there  is  an  (l‘fl)-st  element  and  it  is  Center  n;  and. 
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for  completeness,  (Iv)  the  probability  Is  r(0,Nfl)  .  that  a  routing  will  be 
empty.  It  Is  easily  seen  that  the  probability  of  routing  (n  ,  n  ,  .  .  .  in.) 

j.  2  X 

Is  the  product, 

r(0,nj^)  rCn^^.n^)  ,  .  .  r(Uj^,IH-l)  . 

Two  assumptions  are  made: 

ASSUMPTION  (2.3).  For  each  m  €[0,N]  ,  the  set 

(  r(m,n)  [  n  ^[l,Nfl])  is  a  probability  distribution. 

ASSUMPTION  (2.4).  The  set  of  equations , 

N 

(2.5)  e(n)  a  r(0,n)  +  e(m)  r(m,n)  ,  n  6  [1,N]  , 

has  a  unique  solution,. {  o(n)  |  n€[l,N])  ;  and  the  e(n)  are  all  non- 

negative. 

Assumption  (2.3)  Is  necessary  In  view  of  the  Interpretation  of  the 
r(m,n)  .  In  the  presence  of  (2.3),  Assumption  (2.4)  can  be  shown  to  be 
equivalent  to  the  requirement  that  a  routing  terminate  with  probability  one. 

It  is  Intuitively  obvious  (and  can  easily  be  proved)  that  e(n)  Is  the 
expected  value  of  the  number  of  appearances  of  Center  n  on  a  routing. 

3.  Mathematical  Model 

The  Informal  Interpretations  of  Section  2  are  now  given  mathematical 
substance  by  stating  the  conditional  probabilities  that  the  generic  Jobshop¬ 
like  queueing  system  will  be  in  various  states  j  =  (J,,J_,  .  .  .  ,  J„) 

X  2  19 

at  time  t+h  ,  given  that  it  is  in  state  k  s  •  •  •  >  At 

time  t  .  These  probabilities  are  as  follows: 

N  N  .  « 

1  -  h  \(S(k))  £  ,  r(0,n)  -  h  E  ^J^(n,k  )  (l-r(n,n))  +  o(h)  ,  if  J  =  k  ; 
n=l  n=l  n 

h  X(S(k))  r(0,x)  +  o(h)  ,  if  J  =  k  except  j  =  k  ■!•  1  ; 

X  X 

h  p,(x,k  )  r(x,Nfl)  +  o(h)  ,  If  j  =  k  except  J  =  k  -  1  ; 

X  XX 

b  ia.(x,k^)  r(x,y)  +  o(h)  ,  if  j  =  k  except  *  *^x  ~  ^  *  ^y  *  ^  ’ 

o(h®  ^)  ,  for  J  with  |s(j)  -  S(k)  |  t*  s  >  1  ; 
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where  x  iii  the  second  and  third  expressions  runs  through  [1,N]  ;  and 

X  and  y  In  the  fourth  expression  run  through  [l,N]  ,  omitting  pairs 
with  X  =  y. 

The  first  four  of  the  above  expressions  reflect,  respectively,  the 

possibility  of  no  "events"  occurlng  in  the  interval  between  times  t  and 

t-fh  (except  perhaps  the  completion  of  a  service  for  a  customer  who  then 

requires  another  service  at  the  same  center),  a  customer  arrival,  the 

completion  of  the  last  service  on  a  customer's  routing,  and  the  completion 

Of  one  service  for  a  customer  who  then  requires  another  service  at  a  different 

center.  The  fifth  expression  reflects  the  need  for  at  least  |s(j)  -  S(k)  | 

events  to  intervene  between  the  appearances  of  states  j  and  k  . 

The  process  defined  by  the  proccdint;  expressions,  with  the  notations 

and  subject  to  the  assumptions  of  Section  2,  will  be  referred  to  as 

Jobshop-like  queueing  system  (N,L,M,R),  or  more  briefly  as  System  (H,L,M,R). 

Time-dependent  state  probabilities.  Let  P(k,t)  be  the  probability 

that  System  (N,L,M,R)  is  in  state  k  at  time  t  ;  and  let  P(k,t' |j,t) 

be  the  conditional  probability  of  state  ic  at  time  t'  ,  given  state  j 
at  time  t  .  From  elementary  probability  theory, 

P(k,t+h)  =  2P(j,t)  P(k,t+hjj,t)  , 

the  sum  extending  over  all  state  vectors  j  ;  for  each  state  vector  k  . 

If  the  conditional  probabilities  given  previously  are  substituted  in  this 
equation;  and  then  P(k,t)  is  subtracted  from  both  sides;  and  then  both 
sides  are  divided  by  h  ;  and  then,  finally,  the  limit  is  taken  as  h  tends 
to  zero;  the  final  result  is  the  following  differential-difference  equation: 

(3.1)  =  -[X(S(k))  +  2,*(n,k  )  (1  -  r  (n,n))]  P(k,t) 

QX  u  n 

+  ^  \(S(k)-i)  r(0,n)  P(h(n),t) 

+  ^  l•!'(n,i:^+l)  r(n,N+l)  P(i(n),t) 

^^(n.k^+l)  r(n,m)  P(J(m,n),t) 


+ 


where  the  sums  extend  over  [l,N]  ,  except  that  pairs  with  m  =  n  are 

omitted  from  the  double  sum,  and  terms  are  omitted  for  which  the  vector 

argument  of  P  has  a  negative  component;  h(n)  =  k  ,  except  that  its 

n-th  component  is  -  1  ;  i(n)  =  k  ,  except  that  its  n-th  component 

is  k  +  1  ;  and  J(m,n)  =  k  ,  except  that  its  m-th  component  is  k  -  1 
n  XQ 

and  its  n-th  component  is  k^  +  1. 

4,  Equilibrium  for  Jobshop-Like  Queueing  System  (N,L,M,R) 

By  definition,  an  equilibrium  state  probability  distribution  for 

System  (N,L,M,lO  is  a  probability  distribution,  lp(k))  ,  over  state 

vectors  k  ,  such  that  P(k,t)=p(k)  specifies  a  (constant)  solution 

to  equations  (3.1).  If  such  a  distribution  exists,  it  is  unique,  and 

for  each  state  vector  k  , 

lim  P(k,t)  =  p(k)  , 
t=“ 

e 

the  limits  being  independent  of  the  "initial"  state  of  the  system.  Thus, 
the  equilibrium  state  probability  distribution  provides  a  rather  complete 
description  of  tho  long-run-average  behavior  of  the  state  of  the  system. 
The  theorem  of  this  section  gives  a  condition  under  which  such  a 
distribution  exists,  and  specifies  the  distribution. 


5.  That  the  P(k,t)  approach  limits  follows  from  Theorem  2.9, 
p.  102,  A.T.  Darucha-Reid,  Elements  of  the  Theory  of  Markov 
Processes  and  their  Applications  (McGraw-IIiil,  1900).  It 
can  be  argued,  using  Assumptions  (2.1)  and  (2.2)  that  every 
non-transient  state  of  System  (N,L,M,R)  communicates  with 
every  other,  and  hence  tnat  if  a  limiting  probability 
distribution  exists,  it  is  unique.  The  assertions  then 
follow  from  the  observation  that  the  equilibrium  state 
probability  distribution  is  one  possible  "initial" 
distribution  over  state  vectors. 
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The  following  notations  will  be  used  In  stating  the  theorem,  and 

3 

subsequently: 

K-1 

(4.1)  WdO  =  TT  X(l)  ,  for  K  =  0,  1,  2,  .  .  .  ; 

1=0 


(4.2) 


(4.3) 


w(k) 


T(K) 


IT  k 

=  TT  TT'^[e(n)/tJ.(n,l)]  ,  for  k  =  (k  ,  k  ,  .  . 

n=l  1=1  ^ 

a  state  vector; 

=  2w(k)  ,  svunmed  over  state  vectors  k  with 


for  K  =  0,  1,  2,  .  .  .  ; 


S(k)  =  K 


t 


(4.4)  TT  s  V/(I0  TdOj-  ,  if  the  sum  converges, 

=  0  ,  otherwise. 

It  Is  easy  to  verify  that  the  sum  In  (4.4)  converges  to  a  positive  number 
or  diverges  to  plus  infinity. 

The  proof  of  the  following  theorem  (after  the  preceding  development) 
is  a  straightforward  matter  of  verifying  that  (4.6)  does  define  a 
probability  distribution  ("obvious"),  and  that  equations  (3.1)  are 
satisfied  by  this  distribution  (by  direct  substitution  and  a  moderately 
tedious,  but  routine,  algebraic  reduction): 

THEOSEM  (4.5).  n  >  0,  then  a  unique  equilibrium  state  probability 
distribution  exists  for  Jobshop-llke  queueing  system  (N,L,;<1,R),  and  is 
given  by 

(4.6)  p(k)  =  TTw(k)  V/(S(k))  , 

for  k  a  state  vector. 

An  example.  The  above  theorem  is  extended  in  Section  5,  and  the 
statement  of  remarks  and  derivative  conclusions  will  be  def eri'ed  until  this  is 
done.  As  a  preliminary  illustration  of  the  theorem,  consider  a 


6.  Empty  products  are  assigned  the  conventional  value,  +1 
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two-servlce-center  Jobshop-like  queueing  system  for  which  the  parameters 
are  as  follows:^ 


(4.7) 

X(K)  = 

a/(K  +  1)  ,  v/here 

a  >  0 

and 

X  2^  0  ; 

(4.3) 

Kl,k) 

y 

=  b  k  ,  where  b  > 

0  and 

y  1 

0  ; 

(4.9) 

M-(2,k) 

=  c  k^  ,  where  c  > 

0  and 

z  > 

0  ; 

(4.10) 

r(0,l) 

=  r(l,2)  =  r(2,3) 

=  1  ; 

the 

other  r(m,n)  =  0  . 

The  specification  of  R  implies  that  each  customer  requires  exactly  two 
services,  the  first  at  Center  1  and  the  second  at  Center  2;  and  it  follows 
that  e<l)  =  e(2)  =  1  . 

After  substituting  in  (4.1)  through  (4.4),  it  is  not  difficult  to 
verify  that  tt  >  0  unless  both  (i)  x  =  0  and  (ii)  either  y  =  0  and 
a/b  >1  or  z  =  0  and  a/c  ^  •  If  4?  >  0  ,  then  Theorem  (4.5)  provides 
the  joint  equilibriiun  probability  distribution  of  queue  lengths  at  the  two 
centers,  as  follov's: 


(4.11) 


TT(a/b>^  (a/c)'^ _ _ 

PU,JJ  =  (ijjy 


The  numerical  value  of  tt  can  be  obtained  for  specific  values  of  the 
parameters  by  using  the  fact  that  the  p(l,j)  must  sxim  to  one. 

Variations  upon  this  example  will  be  used  again  in  subsequent 
sections. 


5.  Triggered  Arrivals  and  Service  Deletions 
In  this  section,  the  concept  of  System  (N,L,M,R)  is  generalized 
to  include  cases  where  the  immediate  injection  of  a  new  customer  is 


7.  The  particular  example  used  here  was  suggested  by  R.  \'f,  Conway 

and  V/.  L.  Maxwell,  "A  Queueing  Model  with  State  Dependent  Service 
Rates,"  Joui^nal  of  Industrial  Engineering,  vol.  XII,  no,  2  (1962); 
which  treats  a  queueing  system  equivalent  to  a  single  one  of  the 
example's  two  Interacting  service  centers. 
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triggered  whenever  the  total  number  of  customers  falls  below  a  specified 
limit (  or  where  a  service  Is  deleted  If  a  queue  grows  beyond  a  specified 
maximum  length.  The  generalized  model ,  to  be  referred  to  subsequently  as 
Jobshop-llke  queueing  system  (N,L,M»a)*,  or  as  System  is 

Identical  with  System  (N,L,M,R),  except  for  the  following: 

(I)  There  Is  some  non-negative  Integer,  K*  ;  such  that  the  initial 
state  k  satisfies  S(k)  >  K*  ;  and  if  the  last  service  on  some  customer's 
routing  Is  completed  when  the  system  Is  In  a  state  k  with  S(k)  m  K*  , 
then  the  arrival  of  a  new  customer  occurs  automatically  at  that  Instant. 

(II)  For  each  Center  n,  there  Is  a  k^*  ,  which  may  be  a  positive 
integer  or  +  «  ; such  that  the  Initial  state,  k  ss  (k^,  k^,  .  .  .  »  hj,)  , 
satisfies  k^  <  k^*  for  r  €  [l,Nl  ;  and  if  a  customer  arrives  at  center  n 
when  the  queue  length  there  is  k^*  ,  then  the  imminent  service  required 
by  that  customer  at  Center  n  is  deleted  from  his  routing,  and  he  proceeds 
from  there  immediately,  as  if  he  bad  been  served.^ 

It  is  implicit  in  the  above  statemeiits  that  K*  <  k  ♦  +  k  *  +  ,  ,  .41-* 

-  1  2  •  •  •■» 

If  K*  >  0  ,  then  the  X.(K)  for  K  <  K*  are  obviously  meaningless,  and. 

It  <  +  ai  ,  then  the  M>(n,k)  for  k  >  k^*  are  obviously  meaningless. 

System  (N,L,M,R)  is  the  special  case  of  System  (N,L,U,R)»'  with 
K*  «  0  and  each  k^*  =  4®  .  To  formalize  System  (N,L,M,R)’»  would  involva 
rewriting  the  conditional  probabilities  (beginning  of  Section  3)  for  cases 
where  S(k)  a  K*  and/or  one  or  more  k^  a  k^*  ,  and  rewriting  equations  (3.1) 
correspondingly.  This  will  not  be  done  here,  because  the  modified  equations 
are  of  insufficient  Interest  to  Justify  the  development  of  (necessarily 


Or,  more  generally,  some  customer  is  emitted  from  the  center; 
with  the  restriction  that  the  choice  of  the  customer  to  be 
emitted  must  not  depend  upon  the  future  routing  or  service 
requirements  of  any  customer. 


0. 
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rather  formidable)  appropriate  notations.  ESut  It  Is  necessary  to 
give  generalizations  of  equations  (4.1)  through  (4.4)  to  prepare  for 
the  generalization  of  Theorem  (4.5): 

(5.1)  W*(I0  =  0  ,  for  K  <  K*  , 

K-1 

=  TT  X(i)  ,  for  K  >  K*  ; 
i=K* 

N  k 

(5.2)  w*(k)  =  TT  TT®  [e(n)/M-(n,i)]  ,  if  k  <  k  *  for  n  €  [l,N]  , 

n=l  i=l  n-  n 

=  0  ,  otherwise: 

(5.3)  T*(K)  ss  I.w*(k)  ,  summed  over  k  with  S(k)  =  K; 

(5.4)  TT*  =  {^0  T*(K)j^  ^  ,  if  the  siun  converges, 

=  0  ,  otherwise 

If  K*  55  0  and  each  of  the  k^*  ss  +“  ,  then  equations  (5.1)  through 

(5.4)  specialize  to  equations  (4.1)  through  (4.4),  upon  removal  of 
the  asterisks. 

The  following  theorem,  which  can  be  proved  by  substitution  in  the 
modified  verslori  of  equations  (3.1),  differs  from  Theorem  (4.5)  only  in 
the  presence  of  the  asterisks: 

THEOREM  (5.5).  tt*  >  0  ,  then  a  unique  equilibrium  state 
probability  distribution  exists  for  Jobshop-like  queueing  system 
(N,L,M,R)*,  and  is  given  by 
(5.6)  p(k)  ss  TT*  w*(k)  W*(S(k))  , 

for  k  a  state  vector. 

Remarks .  It  is  natural  to  conjecture  that  if  tt*  =  0  ,  then  no 
equilibrium  state  probability  distribution  exists  for  System  (N,L,M,R)*, 
but  I  have  been  unable  to  prove  this  except  in  special  cases.  If  either 
(1)  there  is  some  (finite)  K  such  that  ?i.(K)  s  0  for  ZC  >  K 


,  or 
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(11)  all  of  the  k  ♦  are  finite;  then  only  finitely  many  terms  In  the 
n 

sum  of  (5.4)  can  be  positive,  so  It  Is  assured  that  nt'  >  0,  and  the 
theorem  applies. 

The  "discovery"  of  Theorem  (4.5)  resulted  from  making  a  sequence 
of  "guesses"  concerning  more  and  more  general  Jobshop-llke  queueing 
systems,  and  proving  successively  more  general  versions  of  the  theorem. 
The  way  in  which  this  took  place  may  be  suggested  by  the  discussion  of 
certain  special  cases  In  Section  6. 

Some  corollary  results.  Suppose  tnat  System  (N,L,M,B)*  satisfies 
the  condition,  tt*  >  0,  so  that  Theorem  (5.5)  applies.  Write  p(S=K) 
for  the  equilibrium  probability  that  the  total  number  of  customers  In  the 
system  Is  K  ;  l.e.,  that  the  state  k  satisfies  S(k)  s  K  .  If 
p(S  =  K)  >  0  ,  and  k  is  a  state  vector,  write  p(k|K)  for  the 
equilibrium  conditional  probability  of  state  k  ,  given  that  the  total 
number  of  customers  In  the  system  Is  K  ;  l.e.,  p(k|K)s:  p<k)/p(Ss:K) 

If  S(k)  =  IC  ,  but  otherwise  p  (k|K)  s  0  .  It  follows  directly  from  the 
theorem  that: 

(5.7)  p(S  =  K)  =  Tt*  T*(K)  W*(K)  ; 

(5.8)  p(k|s(k))  =  w*(k)/T*(S(k))  ,  if  p(S  =  S(k))  >  0  . 

The  probabilities  (5.7)  are  of  interest  in  themselves,  because  there 
are  Instances  where  the  total  number  of  customers  In  the  system  is  of 
primary  concern.  They  also  lead  to  the  overall  mean  arrival  rate  of 
customers  under  equilibrium,  E(X)  : 

(5.9)  E(X)  X(K)  p(S  =  K) 

=  n*  W*(K+1)  T*(K)  . 
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Equation  (5.9)  Is  potentially  of  special  Interest,  because  It  allows 
the  determination  of  "loss  rates"  when  the  variation  in  the  X(K)  Is 
due  to  the  failure  of  "potential  customers"  to  enter  the  system  for 
service  when  too  many  customers  are  already  present. 

The  main  Interest  of  the  conditional  probabilities  (5.&)  perhaps 
stems  from  the  observation  that  they  do  not  depend  upon  the  customer-arrival 
process  —  except  for  its  role  in  determining  whether  n*  >  o  ,  and  poss-lbly 
in  determining  the  values  of  IC  for  which  p(S  =  K)  >  0  .  In  a  sense,  the 
effect  of  particular  customer-arrival-process  parameters  is  concentrated  upon 
the  determination  of  the  p(S  s  K)  . 

Note  also  that  the  parameters  of  the  routing-generation  process 
influence  the  equilibrium  state  probability  distribution  only  by  way  of 
the  e(n)  .  It  was  remarked  previously  that  e(n)  is  the  expected  number 
of  appearances  of  Center  n  on  a  routing. 

Example.  Consider  a  system  identical  with  the  example  given  at  the 
end  of  Section  4,  except  that  for  some  K*,  k,*,  and  k  *  :  (i)  if  the 

total  number  of  customers  in  the  system  falls  below  K*,  then  a  new 
customer  is  instantly  "created";  (ii)  if  a  (k^^*  +  l)-st  customer  joins  the 
queue  at  Center  1,  he  immediately  moves  on  to  Center  2,  without  being 
served  at  Center  1;  and  (ill)  if  a  (k  *  +  l)-st  customer  joins  the  queue 
at  Center  2,  he  immediately  leaves  the  system  without  being  served  there. 
Equation  (4.11)  is  still  valid  for  the  p(i,J)  such  that  i  +  j  >  K*  , 
i  <  kj^*  ,  and  j  ;  but  all  other  p(i,J)  =  0  .  There  is  a  change, 

of  course,  in  the  numerical  value  of  rr  (which  might  be  replaced  by  n* 
for  consistency  with  the  notation  of  the  present  section) . 
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6,  Special  Cases 

The  first  part  of  this  section  is  restricted  to  the  special  case 
of  System  (N,L,M,R)*  in  which  \(K)  =  X(0)  ,  to  be  referred  to  as  the 
constant-arrival-rate  case  of  System  (N,L,M,R)*.  For  this  case,  it  is 

convenient  to  define  the  following  notations,  for  n6  [l,N]  : 

k 

(6.1)  w  (k)  =11  IX(0)  e(n)/ii(n,i)]  ,  for  k  <  k  ♦  , 

n  lajL  —  n 

=0  ,  for  k  >  k  ♦  ; 

’  n 

(6.2)  p  (k)  =  w  (k)/2  w  (i)  ,  if  the  sum  converges, 

“  i=0 

3  0  ,  Otherwise,  for  k  =  0,  1,  2,  .  .  .  . 

The  following  theorem  is  the  direct  specialization  of  Theorem  (5.5)  to 
the  constant-arrival-rate  case: 

THEOREM  (6.3).  If,  in  the  constant-arrival-rate  case  of  System 
(N,L,M,R)*,  Pjj(O)  >  0  for  n6[i,N3  ;  then  a  unique  equilibrium 
state  probability  distribution  exists  for  the  system  and  Is  given  by 
the  product, 

(6.4)  p(Uj^,  kg,  .  .  .  ,  kj^)  =  Pj^(k^)  P2<V  •  •  •  ’ 

for  (kj^,  kg,  ...  ,  liijj)  a  state  vector. 

Interpretation.  When  Theorem  (6.3)  applies,  it  is  immediate  from 

(6.2)  that  the  fixed  n  form  a  probability  distribution 

over  k=0,l,2,  ....  This  distribution  is,  in  fact,  the  equilibrium 
distribution  of  queue  length  for  the  one-service-center  queueing  system 
where  customer  arrivals  form  a  Poisson  process  with  mean  rate  X(0)  e(n)  , 
and  whose  service-completion  process  is  identical  to  that  of  Center  n 
(which  is  well  known,  and  also  follows  from  the  one-service-center  case 


of  the  theorem) 


But  e(ix)  is  the  mean  number  of  appearances  of  Center  n  on  a 
routing  In  the  generic  system  to  which  Theorem  (6.3)  refers,  whence 
X(0)  e(n)  is  in  fact  the  mean  rate  of  customer  arrivals  at  Center  n 
in  this  system.  Thus,  Theorem  (6.3)  can  be  interpreted  as  stating  that 
for  the  constant-arrival-rate  case  of  System  (N,L,M,R)‘i',  the  equilibrirm 
probability  distributions  of  queue  lengths  at  the  individual  centers  are 
Independent;  and  also  each  of  these  distributions  is  identical  with  that 
for  a  one- service-center  queueing  system  "similar"  to  the  center  concerned. 

Further  specializations.  To  represent  an  "ordinary"  M(n)-channel 
Center  n,  with  exponentially  distributed  holding  times  whose  mean  is 
l/P'^  ,  set  k‘‘(n,k)  =  min(k,M(n))  .  If  all  centers  are  thus  specialized 
in  Theorem  (6.3),  the  result  is  the  theorem  of  the  paper  cited  in  footnote  3 
In  this  Instance,  the  sum  of  the  can  easily  be  expressed  in  closed 

form.  Under  this  specialization,  the  one-servlce-ceuter  case  of 
Theorem  (6.3)  gives  the  well-known  equilibrium  distribution  of  queue  length 
for  an  "ordinary"  multi-channel  service  center  with  Poisson  customer- 
arrival  process  and  exponentially  distributed  bolding  times. 

Example.  Suppose  that  in  the  example  at  the  end  of  Section  4,  x  =  0  ; 
and,  to  assure  the  existence  of  equilibrium,  assume  it  is  false  that 
either  y  a  0  and  a/b  >  1  or  z  s  0  and  a/c  2:  ^  *  Then 
P(i,J)  =  Pj^(i)  PgU)  ; 

where  p^d)  =  Pj^(O)  (a/b)V(il)*  and  PgU)  “  (a/c)'^/(ji)^  ; 

the  p  (0)  being  determined  by  the  conditions  that  d.  p  (1)  a  1  , 
n  in 

for  n  =  1,  2  . 

Varying  arrival  rate,  uniform  service  rates.  To  focus  on  the 
effect  of  varying  arrival  rates,  consider  the  special  case  of 
System  (N,L,M,R)*  in  which  e(n)/M'(n,i)  s  1  for  all  1  >  0  .  Then 
w*(k)  a  1  for  every  state  vector  k  .  It  is  easy  to  show  that 
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T*(K)  =  -  1  +  ,  for  K  a  0,  1,  2 . 

It  then  follows  that  a  sufficient  (but  not  necessary)  condition  that 

TT*  >  0  is  that  lim  sup  X(K)  <  1  .  Suppose  that  TT>f'  >  0  .  If  k 

and  j  are  two  states,  with  S<k)  =  K  >  J  =  S(j)  ,  where  p(S  =  K) 

and  p(S  =  J)  are  positive;  then,  from  Theorem  (5.5): 

K-1 

p(k)/p(J)  a  X(i)  ;; 

iaj 

P(S  a  K)  (i  •»  lO  ?v(i)  . 

p(S  =  J)  "  i  +  1 

p(j|K,  =  c  -  i  t  T  ■ 

The  interested  reader  can  improve  his  "feel"  for  the  effect  of  varying 
arrival  rates  by  working  out  more  specialized  examples,  Including  ones 
in  which  p(S  a  X)  >  0  for  only  finitely  many  different  values  of  K  . 

Total  number  of  customers  held  fixed.  Consider,  finally,  the  special 
case  of  System  (N,L,M,R)*  in  which  X(K)  a  0  for  K  ^  K*  .  In  this 
system,  customer  arrivals  occur  when  and  only  when  other  customers  leave 
the  system;  so  the  total  number  present  will  eventually  remain  fixed 
at  K*.  It  follows  from  Theorem  (5.5)  that  p(k)  a  w*(k)/T*(S(k))  if 
S(k)  a  K*  ,  and  otherwise  p(k)  a  0  . 

In  this  instance,  the  model  can  also  be  interpreted  as  representing 
a  closed  system  through  which  the  same  customers  circulate  eternally.  In 
this  interpretation,  the  probability  that  a  customer  served  at  Center  m 
will  next  require  service  at  Center  n  is  T(m,n)  +  r(m,N+l)  r(0,n)  , 
in  the  notation  of  the  present  paper.  If  the  "givens"  of  such  a  system 
include  probabilities  s(m,n)  that  a  customer  served  at  Center  m  will 
next  require  service  at  Center  n,  for  m  and  n  C[l,N];  these 


probabilities  can  be  appropriately  converted  for  the  application  of 
Theorem  (5.5)  by  setting  r(m,n)  =  s(m,n)  for  [l,N-l]  and 
n  €  [1,N]  ,  r(N,N+l)  =  1  ,  and  r(0,n)  =  s(N,n)  for  n  €  [l,N]  . 

As  a  numerical  example,  consider  a  closed  system  made  up  of  one 
"machine  operator"  and  two  "repairmen,"  in  which  there  are  two 
"machines,"  The  machine  operator  functions  like  a  service  center  of 
the  kind  treated  in  the  present  paper;  but  the  "service"  he  renders 
is  to  cause  the  machine  to  "fail."  The  time  to  failure  is  exponentially 
distributed  with  mean  10.  V/hen  a  machine  fails  it  must  be  worked  on 
first  by  one  repairman,  then  by  the  other,  in  a  specified  order.  Each 
repairman  functions  like  a  service  center,  with  mean  service  rate  1/3 
when  he  has  one  machine  to  work  on,  but  with  mean  service  rate  1/2  when 
he  has  two.  Designating  the  machine  operator  as  Center  1,  and  the 
repairmen  as  Centers  2  and  3,  what  has  been  described  is  an  instance 
of  System  (N,L,M,a)*,  with:  N  »  3  ;  K*  =  2  and  \(k)  »»  0  for  K  >  2  ; 
p.(l,l)  a  ^(1,2)  =  1/10  ,  P-(2,l)  =  y-(3,l)  a  1/3  ,  and  Ji-(2,2)  «  M<(3,2)  a  1/2 
and  r(0,l)  =  r(l,2)  =  r(2,3)  =  r{3,4)  =  1  ,  whence  e(l)  =  e(2)  =  e(2>  a  1 
Using  (5,2),  (5.3),  and  the  fact  that  p(k)  =  w*(k)/T*(2)  for  S(k)  =  2  ; 
it  turns  out  that  the  equilibrium  state  probabilities  are  as  follows: 
p(0,0,2)  =  p(0,2,0)  a  6/181  ;  p(0,l,l)  =  9/181  ;  p(l,0,l)  =  p(l,l,0)  = 
30/181  ;  p(2,0,0)  a  100/181  .  It  can  be  concluded,  for  instance,  that 
the  long-run-average  proportion  of  time  during  which  the  machine  operator 
will  actually  have  a  machine  to  operate  is  160/181  ,  or  about  88.4  percent 
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Columbia  University 

Dept,  of  Mathematical  Statistics 

New  York  27,  New  York 

Attn:  Prof.  T.  W.  Anderson 


Applied  Mathematics  Laboratory 
David  Taylor  Model  Basin 
Washington  7,  D.  C. 


Director  of  Operations  Research 
Ernst  &  Ernst 

1356  Union  Commerce  Building 
Cleveland  14,  Ohio 

Attn:  Dr.  E.  Leonard  Arnoff 


Logistics  Department 
The  RAND  Corporation 
1700  Main  Street 
Santa  Monica,  California 

Dr.  Max  Astrachan 


Cornell  University 
Dept,  of  Plant  Breeding 
Biometrics  Unit 
Ithaca,  New  York 

Attn:  Walter  T.  Bederer 


University  of  Washington 
Department  of  Mathematics 
Seattle  5,  Washington 

Attn:  Prof.  Z.  W.  Birnbaum 


Department  of  Statistics 
University  of  California 
Berkeley  4,  California 

Attn:  Prof.  David  H.  Blackwell 


American  Bosch  Division 
American  Bosch  Arms  Corporation 
3664  Main  Street 
Springfield,  Massachusetts 

Attn:  Mr.  R.  K.  Blakeslee 


Armour  Research  Foundation 
Technology  Center 
Chicago  16,  Illinois 


Stanford  University  Dr.  Frederick  Bock 

Department  of  Economics 
Stanford,  California 

Department  of  Engineering 

Attn:  Professor  K.  J.  Arrow  University  of  California 

Los  Angeles  24,  California 


Attn:  Dean  L.  M.  K.  Boelter 
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University  of  Wisconsin 
Department  of  Statistics 
Madison,  Wisconsin 

Attn:  Prof.  G.  E.  P.  Box 


Havard  University 
Department  of  Statistics 
Cambridge,  Massachusetts 

Attn:  Prof.  W.  G.  Cochran 


The  Florida  State  University 
Department  of  Statistics 
Tallahassee,  Florida 

Attn:  Dr.  Ralph  A.  Bradley 


University  of  Maryland 
Mathematics  Department 
College  Park,  Maryland 

Attn:  Prof.  L.  W.  Cohen 


The  George  Washington  University 
Dept,  of  Statistics 
Washington,  D.  C. 

Attn:  Prof.  Harold  Bright 


Headquarters  USAF 
Washington  25,  D.  C. 

Attn:  Mr.  Edward  G.  Bunaway 
Chief,  AFASC-6F 


Bureau  of  Supplies  and  Accounts 
Code  OW 

Department  of  the  Navy 
Washington  25,  D.  C. 


Commanding  Officer 

U.  S.  Naval  Ammunition  Depot  Earle 

Red  Bank,  New  Jersey 

Department  of  Industrial  and 
Engineering  Administration 
Cornell  University 
Ithaca,  New  York 

Attn:  Professor  R.  W.  Conway 


Dunlap  and  Associates,  Inc. 

429  Atlantic  Street 
Stamford,  Connetlcut 

Attn:  Mr.  Gershon  Cooper,  V.  P. 


Case  Institute  of  Technology 
Cleveland  6,  Ohio 

Attn:  Professor  Glen  Cetmp 


University  of  California 
Institute  of  Engineering  Research 
Berkeley  4,  California 

Attn:  Prof.  G.  B.  Dantzlg 


Northwestern  University 
Department  of  Mathematics 
Evanston,  Illinois 

Attn:  Dr.  A.  Chames 


ColiuDbla  University 

Department  of  Industrial  Engineering 
New  York  27,  New  York 

Attn:  Prof.  Cyrus  Derman 


Chief,  Bureau  of  Ships 
Code  732  (A.  E.  Smith) 
Navy  Department 
Washington  25,  0.  C. 


Cornell  Aeronautical  Laboratory,  Inc. 
4455  Genessee  Street 
Buffalo  21,  New  York 


Attn:  Joseph  P.  Desmond,  Librarian 
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General  Electric  Company 
Production  Control  Services 
Room  2401,  570  Lexington  Avenue 
New  York  22,  New  York 

Attn:  Mr.  H.  F.  Dickie 


Director 

National  Security  Agency 
Fort  George  G.  Meade,  Maryland 

Attn:  C3141  (Room  2C087) 


Director 

Operations  Evaluation  Group 
Office  of  Chief  of  Naval  Operations, 
OP03EG 

Navy  Department 
Washington  25,  D.  C. 

CEIR,  Incorporated 

1200  Jefferson  Davis  Highway 

Arlington  2,  Virginia 

Attn:  Dr.  William  Dorfman 


Management  Engineering  Staff 
Bureau  of  Yards  and  Docks 
Navy  Department 
Washington  25,  D.  C. 

Attn:  Mr.  Louis  Doros 


Department  of  Economics 

The  Pennsylvania  State  University 

State  College,  Pennsylvania 

Attn:  Dr.  Lawrence  E.  Fouraker 


jBrown  University 
Department  of  Mathematics 
Providence,  Rhode  Island 

Attn:  Prof.  David  Gale 


General  Electric  Company 
FPD  Technical  Information  Center 
Building  100 
Cincinnati  15,  Ohio 


Section  12.05 

National  Bureau  of  Standards 
Washington  25,  D.  C. 

Attn:  Mr.  Ezra  Glaser 


International  Business  Machines  Corp. 
P.  O.  Box  218,  Lamb  Estate 
Yorktown  Heights,  New  York 

Attn:  Or.  Ralph  Gomory, 

Research  Center 


Cornell  University 

Industrial  &  Engineering  Administration 
Ithaca,  New  York 

Attn:  Prof.  Henry  P.  Goode 


Institute  of  Mathematical  Statistics 
University  of  Copenhagen 
Copenhagen,  Denmark 

Attn:  Prof.  Anders  Hald 


Arthur  D.  Little,  Inc. 

30  Memorial  Drive 
Cambridge  42,  Massachusetts 


Computation  Center 

Stanford  University  Attn:  Mr.  J.  C.  Hetrick 

Stanford,  California 
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Attn:  Prof.  George  E.  Forsythe 


5 


I 


New  York  University 

Institute  of  Mathematical  Sciences 

New  York  3,  New  York 

Attn:  Prof.  W.  M.  Hirsch 


General  Electric  Company 
Management  Consultation  Services 
570  Lexington  Avenue 
New  York  22,  New  York 

Attn:  Alan  J.  Hoffman 


School  of  Industrial  Administration 
Carnegie  Institute  of  Technology 
Pittsburgh  13,  Pennsylvania 

Attn:  Professor  C.  C.  Holt 


Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 

Attn:  Dr.  R.  A.  Howard 


University  of  Minnesota 

School  of  Business  Administration 

Minneapolis  14,  Minnesota 

Attn:  Prof.  Leonid  Hurwlcz, 

Department  of  Economics 


Cornell  University 
Department  of  Industrial  and 
Engineering  Administration 
Ithaca,  New  York 

Attn:  Dr.  Donald  L.  Iglehart 


Industrial  Engineering  Department 
Stanford  University 
Stanford,  California 

Attn:  Prof.  W.  Grant  Ireson 


University  of  California 
Management  Sciences  Research  Project 
Los  Angeles  24,  California 

Attn:  Dr.  J.  R.  Jackson 


The  John  Hopkins  University  Library 
Acquisitions  Department 
Baltimore  18,  Maryland 


Stanford  University 
Department  of  Mathematics 
Stanford,  California 

Attn:  Prof.  S.  Karlin 


Hauchly  Associates  Inc. 

Fort  Washington,  Pennsylvania 

Attn:  Mr.  J.  E.  Kelley,  Jr. 


Cowles  Commission  for  Research  in 
Economics 
Yale  University 
New  Haven,  Connecticut 

Attn:  Prof,  T.  C.  Koopmans 


The  Research  Triangle  Institute 
Statistics  Research  Division 
505  West  Chapel  Hill  Street 
Durham,  North  Carolina 

Attn:  Dr.  Malcolm  R.  Leadbetter 


Department  of  Industrial  Engineering 
The  Technological  Institute 
Northwestern  University 
Evanston,  Illinois 

Attn:  Mr.  R.  N.  I.ehrer 

Professor  &  Chairman 


Mason  Laboratory 
Department  of  Mathematics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York 

Dr.  C.  E.  Lemke 


Stanford  University 

Applied  Mathanatlcs  &  Statistics  Lab. 

Stanford,  California 

Attn:  Prof.  Gerlad  J.  Lleberman 


Columbia  University 

Department  of  Industrial  Engineering 

New  York  27,  New  York 

Attn:  Prof.  S.  B.  Llttauer 


Mr.  H.  D.  McLoughlln 

7709  Beland  Avenue 

Los  Angeles  45,  California 


U.  S.  Army  Chemical  Corps 
Biological  Laboratories 
Fort  Detrick 
Frederick,  Maryland 

Attn:  Dr.  Clifford  J.  Maloney 


Logistics  Research  Project 
The  George  Washington  University 
707  22hd  Street*  N.  W. 

Washington  7,  D.  C. 

Attn:  Dr.  W.  H.  Marlow 


C.  H.  Masland  &  Sons 
Carlisle,  Pennsylvania 

Attn:  Assistant  Controller 


U .  S .  Naval  Training  Device  Center 
Port  Washington,  L.  I.,  New  York 

Attn:  Mr.  Joseph  Mehr, 

Assistant  Planning  Officer 


University  of  Chicago 
Statistical  Research  Center 
Chicago,  Illinois 

Attn:  Professor  Paul  Meir 


£.  I.  DuPont  DeNemours  &  Company,  Inc. 
Engineering  Department 
Wilmington  98,  Delaware 

Attn:  Mr.  F.  F.  Mlddleswart 


Code  250A 

San  Francisco  Naval  Shipyard 
San  Francisco,  California 

Attn:  Mr.  Marvin  0.  Miller 


Dr.  Richard  A.  Miller 
4071  West  Seventh  Street 
Fort  Worth  7,  Texas 


Kaiser  Steel  Corporation 
300  Lcdceside  Drive 
Oakland  12,  California 

Attn:  Hr.  Michael  Montalbano 


Princeton  University 

Dept,  of  Economics  and  Sociology 

Princeton,  New  Jersey 

Attn:  Prof.  0.  Morgenstern 


Corp.  for  Economic  Sc  Industrial  Res. 
1200  Jefferson  Davis  Highway 
Arlington  2,  Virginia 


Attn:  Dr.  Jack  Moshman 
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Department  of  Industrial  Engineering 
Washington  University 
St.  Louis,  Missouri 

Attn:  Mr.  G.  Nadler 


Case  Institute  of  Technology 
Systems  Research  Center 
Cleveland,  Ohio 

Attn:  Dr.  Raymond  Nelson, 
Acting  Director 


University  of  California 
Department  of  Statistics 
Berkeley  4,  California 

Attn:  Prof.  J.  Neyman 


Manufacturing  Controls  Hdqts.  Dept. 
V/estinghouse  Electric  Corporation 
3  Gateway  Center  17-E 
Pittsburgh  22,  Pennsylvania 

Attn:  Mr.  P.  D.  O'Donnell, 

Manager 


University  of  Minnesota 
Department  of  Statistics 
Minneapolis,  Minnesota 

Attn:  Prof.  Ingram  Olkln 


University  of  California 
Department  of  Engineering 
Los  Angeles  24,  California 

Attn:  R.  R.  O'Neill 


President 

Naval  War  College 

Newport,  Rhode  Island 


University  of  California 
Department  of  Economics 
Berkeley,  California 

Attn:  Dr.  Roy  Radner 

^rdue  University 
Department  of  Economics 
Lafayette,  Indiana 

Attn:  Dr.  Stanley  Reiter 


Columbia  University 
Department  of  Mathematics 
New  York  27,  New  York 

Attn:  Prof.  H.  Robbins 


Vulnerability  Analysis  Division 
OASD  (I£cL)  Rm.  30831,  Pentagon 
Washington,  O.  C. 

Attn:  Mr.  Joseph  Romm, 

Staff  Director 


Michigan  State  University 
Department  of  Statistics 
East  Lansing,  Michigan 

Attn:  Prof.  Herman  Rubin 


Brown  University 

Division  of  Applied  Mathematics 

Providence  12,  Rhode  Island 

Attn:  Prof.  M.  Rosenblatt 


Department  of  Philosophy 
Michigan  State  University 
East  Lansing,  Michigan 

Attn:  Prof.  Richard  S.  Rudner 


Attn:  Mahan  Library 
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Colonial  Sugar  Refining  Company t  Ltd. 
1  O'Connell  Street 
Sydney,  Australia 

Attn:  Hr.  R.  W.  Rutledge 


Dr.  Melvin  E.  Salveson 

16  Parish  Road 

New  Canaan,  Connecticut 


University  of  Minnesota 
Department  of  Statistics 
Minneapolis,  Minnesota 

Attn:  Prof.  I.  R.  Savage 


University  of  Michigan 
Department  of  Mathematics 
Ann  Arbor,  Michigan 

Attn:  Prof.  L.  J.  Savage 


Sibley  School  of  Mechanical  Engineering 
Department  of  Industrial  &  Engineering 
Administration 
Cornell  University 
Ithaca,  New  York 

Attn:  Prof.  Andrew  Schultz,  Jr. 


University  of  Pittsburgh 
2930  Cathedral  of  Learning 
Pittsburgh  13,  Penney lviu:i.l  a 

Attn:  Dr.  Alexander  Silverman 


Bureau  of  Supplies  &  Accounts 
(Code  W31) 

Navy  Department 
Washington  25,  D.  C. 

Attn:  Hr.  J.  R.  Simpson 


University  of  California 
Division  of  Electrical  Engineering 
Berkeley  4,  California 

Attn:  Dr.  Otto  J.  M.  Smith 


University  of  North  Carolina 
Statistics  Department 
Chapel  Hill,  North  Carolina 

Attn:  Prof.  Walter  L.  Smith 

Applied  Mathematics  &  Statistics  Lab 
Department  of  Statistics 
Stanford  University 
Stanford,  California 

Attn:  Prof.  H.  Solomon 


The  John  Hopkins  University 
Department  of  Mathematical  Statistics 
34th  &  Charles  Streets 
Baltimore  18,  Maryland 

Attn:  Professor  C.  Stein 


Mr.  Thomas  P.  Styslinger 
Senior  Staff  Industrial  Engineer 
525  William  Penn  Place 
Room  1824 

Pittsbtirgh  30,  Pennsylvania 

Mechanical  Development  Department 
Research  Laboratories  Division 
Bendix  Corporation 
10  1/2  Northwestern  Highway 
Southfield  (Detroit),  Michigan 

Attn:  Mr.  C.  B.  Sung 


Superintendent 

U.  S.  Naval  Postgraduate  School 
Monterey,  California 


Attn;  Library 


Engineering  Library 
University  of  California 
Los  Angeles  24,  California 

Attn:  Mr.  Johanne  E.  Tallman 


Rutgers  -  The  State  University 

Statistics  Center 

New  Brunswick,  New  Jersey 

Attn:  Prof.  Martin  B.  Wilk 


The  Institute  of  Management  Sciences 
Box  62S 

Ann  Arbor,  Michigan 

Attn:  Prof.  Robert  M.  Thrall 


Numerical  Analysis  Research 
University  of  California 
405  Hilgard  Avenue 
Los  Angeles  24,  California 

Attn:  Or.  C.  B.  Tompkins, 
Director 


Princeton  University 
Department  of  Mathematics 
Princeton,  New  Jersey 

Attn:  Professor  S.  S.  Wilks 


Department  of  Mathematics 
Lincoln  Hall 
Cornell  University 
Ithaca  1,  New  York 

Attn:  Professor  J.  Wolfowltz 


Carnegie  Institute  of  Technology 
Graduate  School  of  Industrial 
Administration 
Schenley  Park 

Pittsburgh  13,  Pennsylvania 
Attn:  Prof.  Gerald  L.  Thompson 


Princeton  University 
Department  of  Mathematics 
Princeton,  New  Jersey 

Attn:  Prof.  A.  W.  Tucker 


Professor  Eduardo  Valenzuela 
P.  0.  Box  6643 
Caracas,  Venezuela 


Transportation  Research  &  Development 
Command 

Fort  Eustls,  Virginia 

Attn:  Mr.  Jules  Vichness 
Commanding  Officer 


